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Dynamical Correlation Length near the Chiral Critical Point 
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Abstract. The dynamical evolution of small systems undergoing a chiral symmetry breaking transition 
in the course of rapid expansion is discussed. The time evolution of the dynamical correlation length for 
trajectories passing through a second-order critical point is extracted. It is shown that while the maximum 
value of the correlation length is bound from above by dynamical effects, the time interval during which 
it is near its maximum grows steadily with the system size and with decreasing expansion rate. 

PACS. 25.75.Nq Quark deconfinement, quark-gluon plasma production, and phase transitions - 05.70.Jk 
Critical point phenomena - 64.60.Ht Dynamic critical phenomena - 11.30.Rd Chiral symmetries 



1 Introduction 

The main goal of coUiding heavy ions at high energies is to 
produce matter at high temperature and baryon density. 
In matter at such extreme conditions chiral symmetry may 
be approximately restored (for a review of signatures for 
hot and dense matter see e.g. P). In particular, it has 
been argued that a line of first-order phase transitions in 
the phase diagram of QCD could end in a second-order 
critical point, where the correlation length diverges |2j. 
The critical point in the temperature vs. baryon-chemical 
potential plane has been located on the lattice (HI. Here, 
we study the dynamics of the chiral fields near a critical 
point and determine the behavior of the correlation length 
in finite and rapidly expanding systems, such as the ones 
encountered in heavy-ion collisions. 



2 Chiral Hydrodynamics 

In Chiral Hydrodynamics it is assumed that the long- 
wavelength (classical) modes of the chiral fields evolve in 
the effective potential generated by the thermalized de- 
grees of freedom, which are the matter fields (and possi- 
bly also hard modes of the chiral field). In our model, the 
latter are described as a perfect relativistic fluid, whose 
equation of state is in turn determined by the chiral field 
(via the effective mass), and which can exchange energy 
and momentum with the chiral fields. The chiral symme- 
try breaking dynamics is described by an effective field 
theory, in our case the SU(2)xSU(2) linear cr-model: 



£^q [ij'^df, - g{(T + 75X • tt)] q 
1 

+ 2 



The potential U{<7, tt), which exhibits both spontaneously 
and explicitly broken symmetry, is given by 



[/(a,7r) = ^(a2+7r2-«2)2 



(2) 



Here q is the constituent-quark field q — {u, d). The scalar 
field a and the pseudoscalar field tt together form a chiral 
field (j) — (cr, tt). The vacuum expectation values of the 
condensates are (cr) = /,r and (tt) = 0, where = 93 MeV 
is the pion decay constant. 

For g > 0, the finite-temperature one-loop effective 
potential includes a contribution from the quarks, and is 
given by 



effl 



(3) 



{d^^ad^a + ^^7v^^'^v) - U{a, tt) 



(1) 



Vcs depends on the order parameter field through the ef- 
fective mass of the quarks, ruq = g\(j)\, which enters the 
expression for the energy E. 

For sufficiently small g one finds a smooth transition 
between the two phases. For large coupling, however, the 
effective potential exhibits a first-order phase transition [S] . 
Along the line of first-order transitions the effective po- 
tential has two minima. At T = Tc these two minima are 
degenerate and are seperated by a barrier. As one lowers 
the value of g the barrier gets smaller and the two min- 
ima approach each other. At gc c^: 3.7, finally, the barrier 
vanishes, and so does the latent heat. Below, we study the 
hydrodynamic expansion near this chiral critical point. 

The classical equations of motion for the chiral fields 

are 

d,d^cl>+^-^^0 (4) 

The dynamical evolution of the thermalized degrees of 
freedom (fluid of quarks) is determined by the local con- 
servation laws for energy and momentum. Note that we do 
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not assume that the chiral fields are in equilibrium with 
the heat bath of quarks. Hence, the fluid pressure p de- 
pends explicitly on see for more details. Due to the 
interaction between fluid and fleld the total energy and 
momentum is the conserved quantity: 



denotes the initial size of the system. For the fluid velocity 
we assume a linear profile 



Vr{t = 0,r) = -0{R-r) 
It 



(9) 



(^fluid 



0. 



(5) 



We emphasize that we employ eq. Q not only to propa- 
gate the mean field through the transition but fluctuations 
as well. The initial condition includes some generic "pri- 
mordial" spectrum of fluctuations (see below) which then 
evolve in the effective potential generated by the matter 
fields. Near the critical point, those fluctuations have small 
effective mass and "spread out" to probe the flat effective 
potential. 



3 Correlation Length 

The correlation length provides a measure for the typical 
length scale over which fluctuations of the fields are corre- 
lated. Following the Ginzburg-Landau theory, for an infi- 
nite system in global thermal equilibrium one can expand 
the free energy F arround the thermodynamic expecta- 
tion value 00 - which is given by the minimum of the free 
energy - in a power series: 

F{4>) = ao -f ai4> + a20^ -I- 03^^ + 040" + . • . + an^^ (6) 



Note that the expansion rate (Hubble constant) equals 
l/R, hence larger systems also correspond to less rapid 
expansion. 

The initial conditions for the chiral fields are 



a{t = 0, a;) = Sa{x) + U + 

(-/tt + CToq) • 

Tv{t = 0, a;) = 6tv{x) , 



1 -|- exp 



r- R 



(10) 



a = 0.3 fm is the surface thickness of this Woods-Saxon 
like distribution. Here a^q « is the value of the a field 
corresponding to Ceq. Thus, the chiral condensate nearly 
vanishes at the center, where the energy density of the 
quarks is large, and then quickly interpolates to f^^ where 
the matter density is low. 

(5(7 (a;) and 5tt{x) represent Gaussian random fluctua- 
tions of the fields. We correlate them over an initial cor- 
relation length of ^0 ~ 1-2 fm as described in [Qj. Our fo- 
cus now is on how those "primordial" fluctuations evolve 
through the transition and how the correlation length de- 
pends on the system size. 



where 



and ttN > 0. The correlation length ^ is 4.2 Extraction of Correlation Length 



defined as the second derivative of the free energy 



1 

e 



<fF{4>) 



= 2a2- 



(7) 



4=0 



The correlation length is finite for T since the effec- 
tive potential has a finite curvature there. At Tc the poten- 
tial about 4>Q becomes fiat and thus the correlation length 
diverges as the system approaches the critical point. Of 
course, this can only be true for an infinite system in global 
thermodynamic equilibrium. In high-energy heavy ion col- 
lisions instead one deals with finite and rapidly expanding 
systems where the true dynamical correlation length is 
finite. 



To analyse the time evolution of ^ in our flnite and ex- 
panding system we need to know the dynamical effective 
potential probed by the fluctuating fields. This is done by 
extracting a histogram of the field distribution at every 
time step, within a sphere of radius 1 fm around r — Q. 
This histogram has been averaged over a few random ini- 
tial field configurations, picked according to eq. H10|l . 

The probability distribution is related to the 4-d effec- 
tive action by 



P[cj)] cx exp{-5cff} , 



(11) 



and thus we can extract the effective action " seen" by the 
chiral fields by fitting a polynomial of the form JHl to 



4 Results 



5eff(</')cX-l0g{P[(/.]} 



(12) 



4.1 Initial Conditions 

To study the effects of a finite and expanding system we 
choose the following initial conditions. The distribution of 
the fluid energy density at t = is taken to be spherical: 

e(t = 0,r) = eoqe(i?-r) , (8) 

where Ccq denotes the equilibrium value of the energy den- 
sity taken at a temperature of Tj « 160 MeV which is 
well above the transition temperature Tc ~ 138 MeV. R 



4.3 Numerical Results 

Fig. n depicts the time evolution of the extracted corre- 
lation length ^ for different system sizes. Evidently, the 
correlation length is always finite. For a system of initial 
radius R ~ \ fm, S,{t) ~ ^0 remains approximately con- 
stant during the expansion, equal to the initial correlation 
length. However, for larger R, it develops a maximum at 
intermediate times, roughly twice ^o- Its maximum value 
appears to grow only very slowly with R (once R > 3 fm) , 
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i.e. the approach towards ^ = 00 for i? — ^ 00 is slow. Both 
observations are in hne with who studied finite-time ef- 
fects for infinitely large systems, and estimated that ^ can 
not exceed 2 — 3 times ^q, and that the maximum value of 
^ grows only slowly with decreasing cooling rate. This im- 
plies that experimental observables will not exhibit "fully 
developped" critical behavior even for trajectories through 
the critical point. On the other hand, one might hope for 
some signals to show up even for trajectories that miss the 
critical point somewhat, since the growth of ^ is damped 
by dynamical effects anyways. We also observe from our 
real-time analysis that the time interval during which ^ 
is near its maximum grows steadily with R. This might 
be relevant for observables that integrate over the entire 
collision history. 

The behavior of other couplings at {i > 2) belonging 
to marginal or irrelevant operators will be reported else- 
where. 
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Fig. 1. Time evolution of the correlation length ^ for finite 
and expanding systems of different initial radius R. 



